Correlations between two particles in jets 
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Abstract: We study the correlation between two soft particles in QCD jets. We 
extend the Fong-Webber analysis to the region away from the hump in the single 
inclusive energy spectrum and show that the correlation function should flatten off 
and then decrease for large values of i = In (l/x). 



Introduction 

Perturbative QCD (pQCD) successfully predicted inclusive particle spectra in 
jets (for review see and references therein). It sufficed to make one step beyond 
the leading double log approximation (DLA) and to analyse parton cascades with 
the next-to-leading accuracy by taking into account angular ordering in soft gluon 
multiplication, quark-gluon transitions, exact DGLAP parton splitting functions 
and running coupling effects. The corresponding MLLA (modified leading log 
approximation) "evolution Hamiltonian" has the accuracy 70 + 75% where 70 is the 
characteristic parameter of the pQCD expansion (the DLA multiplicity anomalous 
dimension) given by 
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Here Y = m(Q/2A), with Q the hardness of the jet production process (e + e~ 
c.m.s. annihilation energy) and A the QCD scale. For the LEP-I energy, Q = 91.2 
GeV (Z° peak) we have Y ~ 5.2 (using A ~ 0.25 GeV). 

In spite of expansion parameter 
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being numerically large, 
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the shape of the inclusive en- 
ergy distribution of soft hadrons 
turned out to be well described 
by the analytic QCD curve — the 
so-called limiting spectrum 
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that one derives pushing the min- 
imal transverse momentum in 
partonic cascades Qq down to A. 



Fig. 1 : Inclusive spectrum for Q 

Overall normalization in Fig. 1 (number of hadrons per gluon) is a non-pertur- 
bative parameter that one determines phenomenologically. This parameter should 
cancel in the ratio 
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so that one could expect the two-particle correlation function © to provide a more 
stringent test of parton dynamics. However, for a long time, particle correlations 
are known to be poorly treated by pQCD. 

Here we revisit two-particle correlations and report the results of an improved 
pQCD analysis of the problem l2l . 

Fong-Webber approximation 

The first (and only) pQCD analysis of two-particle correlations in jets beyond 
the DLA was performed by Fong and Webber in 1990. In [3] the next- to-leading 
O (70) correction, C = 1 + ^fal + • • • , to the normalized two-particle correlator 
was calculated. For a system of two quark jets produced in e + e~ annihilation, 

R(h,£ 2 ,Y) = i + i<7 Q (4,£ 2 ,y), 



the Fong-Webber result reads (n j = 3) 
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The first two terms in © are of the DLA origin while the third one constitutes 
an O (70) MLLA correction. This expression was derived in the region \£\ — 
^2\/Y -C 1, that is when the energies of the registered particles are relatively close 
to each other (and to the maximum of the inclusive distribution, see Fig.l). In this 
approximation the correlation function is quadratic in the difference £\ — £2 and 
increases linearly with the sum, £\ + £2- 

Particle correlations from MLLA evolution equations 

To analyse two-particle correlations we write down the evolution equations for 
G Q (*ii^2,y) that follow from the general MLLA evolution equations for jet 
generating functionals EJ. Unlike the case of the inclusive spectrum, equations for 
D^> are inhomogeneous due to the presence of the product of one-particle distri- 
butions D{£i)D{£2). In the small- x limit, by approximating the energy fraction 
integrals of the parton splitting functions as follows, 

d£' P{x') F{£ - H!) ~ / d£' [ a - c 2 5(£') ] F(£ - (!) , 
'0 Jo 

it is straightforward to reduce the original integral equations to a system of two 

linear (inhomogeneous) second order differential equations in variables £\ and 

2/2, where y = ln(A;^/A) and we keep the ratio of particle energies fixed, £\ — 

£2 =const. 

We then substitute the product C ■ D(£\)D{£2) for and solve the equations 
iteratively, using the fact that the normalized correlator is a slowly changing 
function as compared with the distributions themselves. The range of applica- 
bility of the solution so obtained is not restricted to the vicinity of the hump as 
in [3|. Actually, the new solution can be trusted for arbitrary values of £i, as long 
as £1, £2 > 2 (to respect the adopted soft approximation, x <^ 1). 

For example, the answer for the two-particle correlation inside a gluon jet reads 
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The function A is given by 

A : 
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where tpi = h\D{t L ,Y), and the subscripts £ and y mark its respective derivatives. 
Since logarithmic derivatives of the inclusive spectrum are, typically, of the order 
of ipi ~ ip y = O (70), we have A = O (1). This term is already present in the 
DLA, see 0, and is responsible for the fall-off of the correlation with increase of 
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The term b (ipx/ + tp2,t)> where b 

numerator of @ is the next-to-leading (MLLA) correction O (70). The term 8\ is 
given by 
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It is O (70) and constitutes a MLLA correction as well, since xe ~ Xy 
The coiTection term 62 is given by 

^2 = Jo* (XiXy + Xly) , Xly ~ = C (7o) • 
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It is O (7q) and constitutes a NMLLA correction, as well as the term a\i- Thus, all 
terms in the square brackets in © are of the order of 7g. Being formally negligible 
within the MLLA accuracy, we nevertheless keep them in the answer. By so doing 
we follow the logic that was used in the derivation of the single inclusive MLLA 
spectra namely, solving exactly the approximate (MLLA) evolution equations. 
The correlator Cq itself enters on the r.h.s. of Q via © and (0. Substituting 
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into the expressions for 5i and 82 provides then an iterative solution. The sum of 
the correction terms 5\ + 82 +axi is very small, O (lO -2 ) , near the hump (l\ ~ £2) 
and increases away from it. 

The correlators for the quark jet and for e + e~ annihilation (two quark jets) can 
be obtained from that of the gluon jet. 

Below the results are presented for four different bands of the {£\, £2) plane, in 
comparison with the Fong-Webber approximation and the OPAL data [6]. In the 
left column we show R — 1 as a function of the sum £\ + £2 for three values of the 
difference, and on the right column, vice versa, as a function of the difference for 
fixed sum. 
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Conclusions 

The correlation should be strongest when the two particles have the same en- 
ergy, £\ = £2, in agreement with the Fong-Webber analysis. 

At the same time, the correlation function that we derived directly from the 
evolution equations without expanding the answer in the formal perturbation pa- 
rameter 70, no longer increases linearly but flattens off (and then tends to decrease) 
with £\ + £2. Such a behaviour is in accord with general theoretical expectations. 
Indeed, the smallest energy gluon with £ — > Y is pushed at large angles, ~ 1, 
and, in virtue of the QCD coherence, should be radiated independently of the rest 
of the parton ensemble. 

Though our curves are somewhat closer to the OPAL measurements than the 
Fong-Webber results, the discrepancy remains substantial. Whether this discrep- 
ancy is due to higher order perturbative corrections or rather due to non-trivial 
hadronization effects that have not been seen in the inclusive one-particle spectra, 
remains to be studied. Forthcoming experimental data on two-particle correlations 
in g/q jets produced in pp collisions (CDF) should elucidate this problem. 
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